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ABSTRACT 



Using the decoherence formalism of Gell-Mann and Hartle, a quantum system is 
found which is the equivalent of the classical dissipative chaotic Duffing oscillator. 
The similarities and differences from the classical oscillator are examined; in partic- 
ular, a new concept of quantum maps is introduced, and alterations in the classical 
strange attractor due to the presence of scale-dependent quantum effects are stud- 
ied. Classical quantities such as the Lyapunov exponents and fractal dimension are 
examined, and quantum analogs are suggested. These results are generalized into a 
framework for quantum dissipative chaos, and there is a a brief discussion of other 
work in this area. 
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I. INTRODUCTION 



Since classical chaos first began to be studied, a conspicuous puzzle has been how to 
reconcile this nonlinear, purely classical phenomenon with an underlying linear quantum 
theory. If we believe, as we must to be consistent, that all of physics is fundamentally 
quantum mechanical in nature, then we must further believe that true chaotic systems do 
not exist. At some point, at length scales determined by Planck's constant, the deterministic 
uncertainties of classical chaos must give way to the probabilistic uncertainties of quantum 
mechanics. 

But tackling these problems is not simple. The nonlinear equations of chaos are, in gen- 
eral, only solvable with modern high-speed computers, and their quantum analogs share this 
limitation. Also, chaos itself encompasses two major types of behavior: Hamiltonian chaos 
in which energy is conserved, and dissipative chaos. It is in principle straightforward to find 
and solve quantum equivalents to Hamiltonian systems, if difficult in practice, and consid- 
erable progress has been made in recent years in understanding these systems. Dissipative 
systems are much more foreign to quantum mechanics as it is usually studied. 

Recently, Murray Gell-Mann and James Hartle have used the decoherence functional 
formalism of quantum mechanics to show how quasiclassical laws can arise from an under- 
lying quantum theory I applied this approach to the problem of Brownian motion, 
demonstrating how their scheme reproduces exactly the classical Langevin equation in a 
fairly broad class of systems @]. 

A natural next step is to apply this to systems with interesting classical behavior. Since 
dissipation is easily and indeed naturally included in such systems, an obvious candidate 
for study is dissipative chaos. Once a quantum system is found whose limiting behavior is 
equivalent to a classical chaotic system, we can study how the residual quantum mechanical 
effects alter the system, and what difference this makes to the classical behavior. 

In section II we briefly examine the family of quantum systems from which we will draw 
our model, and derive the quasiclassical equations of motion for them. We then go to the 
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limit of an infinite reservoir of oscillators with a continuum of frequencies, and specialize to 
the case of a forced, damped nonlinear oscillator. 

In section III we examine the classical behavior of one such system, the damped, driven 
Duffing oscillator. There is a brief discussion of dissipative chaos, the structure of the strange 
attractor, and the bifurcations leading to chaos. Several quantities useful for characterizing 
the chaotic behavior are defined: the fractal dimension and Lyapunov exponents, and their 
relationships are examined. 

In section IV we look at the decoherence functional, and define the idea of a quantum 
map. The system is examined as a Wigner distribution, and we see how the invariant measure 
of the strange attractor goes over to the quantum case. Problems of coarse-graining and 
decoherence are discussed. Then we look at the system from a Master equation point of 
view, and compare this description to the decoherence functional approach. In section V 
we see how the various classical quantities used to characterize chaotic behavior can be 
reinterpreted for our quantum system, by treating it as a classical system with noise for 
sufficiently coarse length scales. 

A few other treatments of quantum dissipative chaos are mentioned in section VI, and the 
differences between Hamiltonian and dissipative chaos are pointed out. Finally, in section 
VII a case is made for a general theory of quantum dissipative chaos. 

II. DAMPED DRIVEN QUANTUM SYSTEMS. 

A. The Quantum Systems. 

Picking a good set of candidate systems requires some thought. Many widely studied sets 
of chaotic equations have only a loose connection to actual physical systems; many others 
are extreme coarse grainings of very complicated systems with many degrees of freedom, 
e.g., fluid dynamics. It is much better to deal with comparatively simple systems, whose 
decoherence functionals can be calculated easily. For this reason, I have elected to study 
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damped driven nonlinear oscillators, which can be easily modelled as particles moving in a 
potential well, interacting with a reservoir of simple linear oscillators. In particular, I will 
concentrate on one such system, the damped, driven Duffing oscillator. 

Earlier work has chiefly considered systems with reservoirs in a thermal state. For the 
purpose of this model, I wish to consider instead a system whose reservoir is initially in a 
coherent state. 

Consider a system of harmonic oscillators. We assume them to be in a state 
where {z/} represents a set of complex numbers uj. If cij is the annihilation operator for 
the ith oscillator, then aj|{z/}) = i^jKz/}). 

As shown in earlier papers 0,^ , the decoherence functional for a system interacting with 
a reservoir is 

D[x'{t),x{t)] = exp|z(5sys[2;'(t)] - S^yMt)])/^] J SQ'5Q 5{Q'{tf) - Q{tj)) 

Xexp|z(5res[Q'(t)]-5,es[g(t)] 

- {V{x'{t), Q'{t)) - V{x{t), Qmdt)/n^ 
xp{x'o,Q'q;xo,Qo) 
= exp|z(S'sys[x'(t)] - Ssys[x{t)] + W[x'{t),x{t)])/n^p{x'o;xo). (2.1) 

Here 5'sys[x(t)] is the action of the system for a given trajectory x(t), SrcslQit)] is the action 
of the reservoir for a given trajectory Q{t), and V{x, Q) is the interaction potential between 
the system and reservoir variables. We will assume this to be a bilinear potential of the form 

V{x,Q) = -J2jkxQ'', (2.2) 

k 

where is the coordinate of the kth oscillator in the resevoir. We will shortly allow the 
number of oscillators to go to infinity, and assume a continuum of oscillator frequencies, but 
for now let us deal with the discrete case. 

We will also assume that the density matrix factors: 

pix'o, Q'o, xo, Qo) = xix'o; xo)0((5'o; Qo), (2.3) 
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where (f>{Q'o',Qo) = {Q'o\{^}){{^}\Qo) is the pure coherent state described above. 
We can readily calculate the influence functional for this system. It is just 

exp|ziy[x'(t),x(t)]A} = 1 1 j dQfdQodQ',K:,^,^{Qf,Q',)K,^t^{Qj-Qo)<P{Q'o,Qo) 

= {{^}\Slit)Sm\{^}). (2.4) 

where K^(j^^{Qf] Qq) is the transition amplitude from Qq to Qf of the reservoir and S^^t) is 
the time evolution operator of a forced harmonic oscillator driven by the time-dependant 
interaction V{x{t), Q) given in ( |2.2| ). This is a well-known problem For a single oscillator 
of frequency uj the operator is 



S^^t) = exp 

where 



aa' — a a 



D{a), (2.5) 



y/2mujh J to 

and 



a = ^ r e''^'xis)ds, (2.6) 



D{a)\iy) = \iy + a)e''"'^'"' >. (2.7) 

For 1/ = this just reduces to the usual form of the influence functional for an oscillator 
initially in the ground state: 

expliW[x'{t),x{t)]/?i\ = exp —l + -p- — f dt f ds cos{uj{t-s)){x'{t)-x{t)){x'{s)-x{s)) 



h [ Amu 



k) ''to 



dt J ds sin(cu(t - s)){x'{t) - x{t)){x'{s) + x(s))|. (2.8) 



2muj Jto Jto 

For non-zero u we get an additional exponent of the form 



Im(2i^a* + 2u*a') = -fJ^-r dt (Reu cos (out) + lmz/sin(tut))(x(t) - x'(t)). (2.9) 

V muh Jto 

Generalizing this to many oscillators, we get the influence phase 

W[x'(t),x(t)] =y-^ dt dscos{iUk{t - s))ix'it) - xit)){x'{s) - x(s)) 
V 4:muJk Jto Jto 



2mujk 



f ^ dt f dssm{ujkit - s)){x'{t) - x{t)){x'{s) + x{s)) 

to J to 



+ 7fc\/ / dt(Rei'kCos(uJkt) +lTah'kSm(uJkt))(x'(t) — x(t)). (2-10) 

V mwk Jto 



For practical purposes, we generally assume that the interaction began at to = and con- 
tinued up until some final time i/, so as to avoid having infinite limits in the integrals. 
We will now assume that the action of the system variables is of the usual form 

-^sysfa^W] = L{x{t),x{t))dt. (2.11) 

Jto 

We can then change variables to 

X^hx + x'), (2.12a) 

^ = x-x', (2.12b) 

and write the decoherence functional in terms of the new variables. This is easily shown to 
be 
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dL_ 
dX 



(X(t),X(t)) 




+lk\ (Rei/fcCos(t^fc(t)) + Imz/fcsin(u;fe(t))) 



oL ~] 



/ dt ds cos{uJk{t — s))^{t)^{s) 




ds sm{ujk{t — s))X{s) 



(2.13) 



Note that the real part of the phase includes the Euler-Lagrange equation of motion for the 
system, with the addition of a retarded force due to the interaction with the reservoir. The 
imaginary part is strictly non-negative, with a minimum at ^{t) = 0, and hence tends to 
suppress D[x'{t),x{t)] for large ^. This makes our expansion in terms of C,{t) seem reasonable, 
and also causes the system to decohere, at least approximately, since ^(t) ^ corresponds 
to off-diagonal terms. The term occurs because of an integration by parts. 



We can now look at the classical system equivalent to the above quantum system, i.e., 
with the same action functional and distribution of oscillators. A coherent state is often 
characterized as a more "classical" state of an oscillator than the usual Fock states; it can 
be thought of as the state of an oscillator begun at a given initial position and momentum, 
within the limits imposed by the uncertainty principle. 

We will begin by assuming knowledge of the trajectory of the system variable x(t), and 
ask what the behavior of the reservoir of harmonic oscillators will be . Assume that we 
start the oscillators in a definite state Q^{t = to) = q^,Q^{t = to) = . The interaction 
potential is linear, so we can treat the trajectory of the system variable x{t) as a simple 
driving force, giving us an equation of motion for the kth oscillator 



B. The Classical Equivalent. 




(2.14) 



The solution to this equation is simply 
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Q*^(t) = cos{uJk{t - to)) + (f V^fc) sm(cjfc(t - to)) 

H / sm{ujk{t — s))x{s)ds. (2.15) 

muk Jto 

If the system is described by a Lagrangian L{x,x), then we can write down the Euler- 
Lagrange equation 

rj 8T 8T 

-^(x(t),x(t)) - ^{x{t),x(t))+Y.^,Q\t) = 0. (2.16) 



We can clearly substitute the above expression ( p.l5|) for Q^{t) in the Euler-Lagrange equa- 
tion to get 

^^(a;W,a;(t)) - -^{x{t),x{t)) \ ^-i^,[q^ co^{uJk{t ~ U)) + (t^V^fc) sin(^fc(t - to)) 



+ 



Ik 



sin(cjfc(t - s))x(s)rfsj = 0. (2.17) 



This expression is clearly closely related to the real part of the phase in the decoherence 
functional, if we make the identity 



1h 

Rei/fc = g^ (2.18a) 

Imi/fc = v^. (2.18b) 



'2huJk 



m 



If we write the above classical equation as e{t) = 0, then the real part of the phase is just 

/*' C{t)e{t)dt. 



How do we interpret the imaginary part of the phase, however? In treating reservoirs in 
an initial thermal state, we identified this term as the effect of a stochastic force F{t) arising 
due to thermal noise. However, in this system, there is no noise classically. The persistence 
of this term indicates a fundamental difference between the quantum and classical systems. 
As Gell-Mann and Hartle point out [Q, in the quantum system there is always noise from 
zero-point oscillations, unlike classical oscillators. So the actual equation of motion derived 
from the quantum theory is 
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° ^ m^''^^^'''^^^^ ~ ii(x(t),x(t)) - F{t) 

+ J2'yk((l''cos{uJk{t - to)) + {v^/u^k) sin(cjfc(t - to)) — f sm{uJk{t - s))x{s)ds] 

u V muk J to J 



(2.19) 

where F{t) is a stochastic force with {F{t)) = 0, and a two-time correlation function 

(Fims)) = E ^ cos(.;.(t - .)). (2.20) 



C. The Continuum Limit. 

In order to consider the sorts of classical systems we are concerned with, we must go to 
the limit of a continuum of oscillator frequencies, both classically and quantum mechanically. 
In doing this, we replace our sums over oscillators with integrals over a distribution function 
g{uj). The usual choice for such a g{Lj) is the Debye distribution 

g{uj) = (2.21) 
[O, uj>VL. 

The reservoir degrees of freedom will become a continuum, Q^{t) — >■ Q{uj,t), and the eigen- 
values Uk will become a continuous complex function i>{u!). In general, Q must be taken to 
be fairly large. More precisely, we want Q ^ l/{tf — to), so that the relaxation time of the 
reservoir is much less than the time-scale of the problem. 

Let's consider now the various components of W[X{t),C,{t)] one at a time. In the con- 
tinuum limit, we have 



V f ds sm{ujk{t - s))X{s) ^ — r duj f ds^^ sin(cu(t - s))X{s). (2.22) 

^ mCUk Jtn 1TL JO Jtn UJ 



We can invert the order of integration and do the u integral, substituting (|2.21|) for g{uj): 



1 ft gi^) V /"* 

— ds duj sin(u;(t — s))X(s) = — — I ds I dujusiniuit — s))X{s), 

m Jto Jo u) m\ V Jto Jo 



r] ft ^ f Qcos{Q{t-s)) , sm{Q{t-s)) \^^ ^ 

— — + (t-s)^ r^^^' 

V r^,^(^f!!!W^W(,). (2.23) 



mfi^ Jto ds\ t — s 



Now we use the fact that Q is large. This imphes that the sinusoidal terms will oscillate 
very rapidly, so that the integral will tend to cancel out to zero. We expect the largest 
contribution to come in the region where s is close to t. Thus, we expand X{s) about t to 

get X{s) X{t) — X{t){t — s) -\ . Substituting this into the above integral, we can solve 

it term by term to get 



- ( ds I du;^sm(u(t-s))X(s)^-^X(t)-^^X(t) + 0(n-^), (2.24) 
m J to Jo UJ mil 2miV 

where the additional terms become small in the limit of a large cutoff Q. 



The second term has the form of a dissipation with constant 2r = 7rr]/2mMQ'^. The 
first term is a linear restoring force. If the system were a harmonic oscillator, this would 
cause a shift in the oscillator frequency. We can absorb this term into the system action as 
an additional harmonic oscillator potential: 

rtf 
'to 

If our Lagrangian is the usual L{x,x) — ^Mx^ — U{x), then we effectively have a new 
Lagrangian 

L{x, x) L\x, x) = ]-Mx^ - U'{x), (2.26) 

where 

In subsequent analysis, it will be U'{x) that we are interested in, as the effective potential. 
The imaginary part of ^] is also of interest. Here we have 

r dsco^{uj^{t-smt)i{s) Tdco f'ds^-^sm(u;{t-smms), 

^ ^mujk Jto 4m JO Jto uJ 

77 /■*/ 

= ds du ujcos{u{t - s))^{t)^{s), 
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If we do this derivative we see 

d /cos(fi(t - s)) - 1\ _ nsm{n{t - s)) cos(^](t - s)) - 1 
d~s[ t^s ) ^ ir~s ^ {t - s)2 ' 

where the first term is larger than the second by a factor of roughly Q{t f — to) ■ In the limit 
of large Q this term approaches a delta function. We again expand ^(s), 



and substitute it into ( |2.28| ) to get 



1] f^f J /^'^Os(r2(t — s)) — ^^^^^^ ^ TCT] \ /n^n-2\ 



,^a^J, ''TA \-s )m^i')''^e(t) + 0(a-'). (2.29) 

We define a new constant K = 7Tr]/2'mQ. 

At this point, someone will likely cry foul. In the real part, we kept the two lowest order 
terms, while in the imaginary part we kept only one! I will give three different arguments 
why this is legitimate: 

1. The real and imaginary parts of W[X{t),C,{t)] serve different purposes in the decoher- 
ence functional, and have different effects. A small correction to the imaginary part has only 
a minor effect on the level of decoherence. The inclusion of a small amount of dissipation in 
the real part, by contrast, leads to qualitatively different solutions. 

2. The leading order term of the real part can be absorbed into the system action by 
going to an effective potential U'{x), as we've seen, so that we must go to the next order to 
observe fundamental changes in the classical equation of motion. 

3. The imaginary part of W[X{t),^{t)] corresponds to a small stochastic noise. In the 
above limit the two-time correlation function becomes {F{t)F{s)) = {7irih/A'mQ)6{t — s) = 
hK6{t — s). In the usual limit where h is small, any corrections to the above correlation 
function would be too tiny to matter. The dissipative term lacks this factor of h. 

Finally, we have the terms arising from the initial condition of the reservoir variables. 
When we go to the continuum limit here, the discrete sum in ( p^.lOD becomes an integral. If 
we choose i^{uj) ~ 6{iJ — luq), then we get 
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/ 2h 
k V f^^k 

where q can be set arbitrarily by adjusting the amphtude of Rez/(c(j). This term has the 
form of a periodic driving force. We can set the phase 0o to zero by making Imz/(ct;) = 0. 
Physically, the presence of this term is equivalent to the system being driven by a plane 
wave at the frequency uq. 

Thus, our decoherence functional becomes 



I ( rf ( ■■ dU ■ \ 

D[X(t),e(t)] =exp-|y^ dt[-MX{t) - —{X{t)) - 2MTX{t) + MqcosMj^it) 

-MXo^o + ^K f^' dt eit) + O(e')}x(x'o; x,). (2.31) 



which gives us a quasiclassical equation of motion 

x + ^'^{x) + 2Tx = qcos(ujot) + F(t)/M, (2.32) 
M dx 

where F{t) is a stochastic force with {F(t)) = and {F{t)F{s)) = hK5{t — s). In a 
completely classical derivation, of course, this stochastic force would be absent. Thus, we 
have found a quantum system equivalent (in the appropriate limits) to a classical nonlinear 
oscillator with a periodic driving force and dissipation. All that remains now is to specialize 
to a chaotic example. 

III. THE DAMPED, DRIVEN DUFFING OSCILLATOR 



The quasiclassical equation of motion ( 2.32| ) is a fairly general expression for a one- 



dimensional damped, driven system. Many such systems exist which exhibit chaotic behavior 
for some values of the constants F and q. The ordinary pendulum is an example of such a 
system, where U'{x) = — cos (a;). We will be examining another system: the damped, driven 
Duffing oscillator. This nonlinear oscillator has a polynomial potential 

— U'(x) = -x" - -x\ (3.1) 
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We choose units to set M = 1. This system has the advantage of having been thoroughly 
studied and examined in the past, and also, since the potential is a polynomial, of not having 
an infinite number of nonzero deriviatives. We will later see that this is convenient, though 
not vital. 

The equation of motion is now 

x + 2Tx + (x^ -x) = q cos(tuot) + F{t). (3.2) 

The potential is double-welled (see fig. 1). For some values of the constants, the oscillator 
undergoes periodic motion. By adjusting the frequency, one causes the system to undergo a 
series of bifurcations until eventually it enters into a region of chaotic behavior, typified by 
the presence of a strange attractor (see fig. 2). If one adjusts the driving force further and 
further the chaotic region is eventually left, and periodic motion returns 0. 

It is convenient to look at the long-term chaotic behavior in terms of a constant phase map 
or surface of section. That is, we consider the position x and momentum p at the discrete 
times ti = 2TTi/ujQ. By then plotting the values Xi and Pi we make the fractal structure of 
the strange attractor very clear (see fig. 3). We are also able to bring the mathematical 
toolbox of discrete dynamical system theory to bear on the problem. Formally, we define 
the constant phase map Xi Xi+i = fx{xi,pi), pi pi+i = fp{xi,pi), where f is an operator 
which evolves the point {xi,pi) in phase space forward in time by 2tx /uoq. 

We can now define a probability measure P{x,p) on our phase space. In this discrete 
dynamics, it evolves according to the equation 

Pi+i{x,p) = J dx j dp'5{x - fx{x',p'))5{p- fp{x',p'))Pi{x',p'), (3.3) 

It is very useful then to consider an invariant measure, which gives a probability distribution 
on the strange attractor. This is defined by the equation 

Pi+i(x,p) = Pi{x,p) = Pinv{x,p). (3.4) 

For a chaotic system such as the Duffing oscillator, Pinv will not be an analytic function; 
rather, it will be a generalized function. Also, it will not in general be unique; there are 
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many invariant measures, most corresponding to unstable solutions, fixed points or periodic 
points. It has been shown that the inclusion of a small amount of noise removes both of 
these objections, eliminating the unstable solutions and making the function analytic 0. 
The inclusion of noise effectively broadens the 5-functions in ( |3.3| ), making it impossible for 
probability measure to remain poised on an unstable fixed-point solution or periodic orbit, 
leaving only the strange attractor as a stable set. One can therefore define the classical Pmv 
as the limit of this unique Pmv as the noise goes to zero; it will still be a generalized function 
with a fractal structure, but now unique. Note that in the quantum system, the noise is 
always non-zero; we will see how this modifies our definitions in sections IV and V. 

The structure of the strange attractor arises as a limit of repeated stretching and folding 
of phase space. As we look closer and closer at the component points of the attractor, we 
see repeated layers of substructure at every scale (see fig. 4). Such infinite substructure is 
commonly characterized by its fractal dimension. There are a number of ways of defining 
dimension, each of which has slightly different properties. 

One common definition is that of the capacity or Kolmogorov dimension, Dq- This is 
calculated by means of a box-counting algorithm. Phase space is divided into small cells of 
linear size e, and one counts the number of cells A^(e) which contain points of the attractor. 
The dimension is then 

IniV(e) , , 

Dc = lini 3.5 

Ine 

Though this definition is fairly easy to calculate numerically, it does not reflect the fact 
that an orbit may visit regions of the attractor with varying frequency. To take this into 
account, one may instead use the information dimension, Dj. Again, phase space is divided 
into cells of linear dimension e. The probability that a given point will fall in the ith cell is 
Pi. Dj is then 

T,Pi In Pi 

D/ = lim^- . (3.6) 

e^o Ine 

If Pi is equal for all cells that are visited, then Dj = Dc', otherwise Dj < Dq. 
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If one examines a small cell of phase space evolve according to our equation, it will tend 
to be stretched along one dimension. The overall phase space volume, however, will contract, 
due to the effects of dissipation. This stretching is what provides the well known signature 
of chaotic systems, sensitivity to initial conditions. The contraction, together with the more 
global process of folding, is what leads to the fractal structure of the attractor. One can 
average these two effects over the entire attractor to calculate the Lyapunov exponents. In 
our two-dimensional phase space, this will be a pair of numbers Ai and A2, with Ai > 
characterizing the stretching and A2 < characterizing the contraction. Since overall phase 
space volume is shrunk by this system, clearly Ai + A2 < 0. If our phase space was n 
dimensional, there would clearly be n characteristic exponents. 

Lyapunov exponents are calculated by considering the time evolution of an infinitesimal 
frame of basis vectors in phase space. One can perform a Gram-Schmidt orthogonalization, 
separating out the most rapidly increasing direction from less rapidly increasing directions 
repeatedly until one has n orthogonal vectors. One then takes the logarithm of the rate of 
change in each of those directions. Allowing the frame to evolve for many driving cycles 
lets one follow a phase space cell as it samples all parts of an attractor. In this way one 
calculates the average values of the exponents. The values of the A, are global quantities, 
characterizing the attractor as a whole, or equivalently, the long-term behavior of orbits 
throughout the attractor. 

Calculating the highest exponent, Ai, is not very difficult. Finding values for a full 
spectrum of exponents, however, is rather tricky, and requires a subtle touch. I refer those 
interested to the papers of Wolf, et al. and Brown, et al. for details [p|,p!0|. These definitions 
for fractal dimension and Lyapunov exponents run into trouble in the quantum case, where 
taking limits as e — > is not very well-defined. We will see in section V how one can adjust 
these definitions appropriately. 

IV. QUANTUM MAPS AND STRANGE ATTRACTORS 
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A. Quantum Maps. 



In the classical case we went from continuous to discrete dynamics by going to the 
constant phase map of the Duffing oscillator. In considering a quantum equivalent, it is 
convenient to coarse-grain our selected trajectories x{t) and x'{t) by considering only their 
values at the times ti = 2'Ki/ujQ of constant phase. The decoherence functional then becomes 

D[{x,},{x'.^]= [ 6x [ 5x'D[x{t),x'{t)], (4.1) 



where the decoherence functional on continuous trajectories D[x{t) , x' (t)] is given by ( [2.31 ). 



The path integrals are over all paths of x{t) and x'{t) which pass through the points Xi and 
x'i respectively at times ti. 

Such a coarse-graining is discussed by Gell-Mann and Hartle 0. In general, in order 
for such a system to be sufficiently decoherent, we must also coarse-grain on the positions 
{xi} and {x'i}. Instead of specifying the positions exactly, we instead require just that the 
positions fall in one of a group of short intervals A^. at the times ti. A history is then given 
by specifying the sequence of a^'s. We'll use the short-hand notation a for this sequence. 

We can estimate the minimum length d of such intervals by requiring that the off-diagonal 
terms of the decoherence functional be strongly supressed for \xi — x'i\ = > d. As Gell- 
Mann and Hartle show, this depends on the separation between times ti, the strength of the 
coupling, and so forth. A rough estimate gives 

For the quantum maps the interval is At = 27r/u;o- ^ is the frequency cutoff; it is equal to 
K/MT. So we get an estimate 

d^ . (4.3) 



Since equation ( |2.31[ ) is expressed in terms of the variables X{t) and ^{t), it might be 



useful to change variables in our coarse-grained systems. As the decoherence functional is 
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supressed for large ^, we can we can treat our integrand as being quadratic in ^{t), and carry 
out the ^ integration. This gives us 



(4.4) 



where e{t) = is the classical equation of motion as given above in (|2.32|) and w{Xq,po) is 
the initial Wigner distribution, obtained by the integral over ^o- 

We see that if the {Xi} do not lie along a classical trajectory e{t) = 0, then the functional 
will be supressed. So the most probable histories are those which lie along the classical 
trajectory. The X path integral cannot be done exactly in most cases, but one can see that 
in general the {Xi} must lie near the {xj} for some classical problem for the probability to 
be of reasonable magnitude. 

The Wigner distribution w{X,p) (not to be confused with the influence phase 
W[X(t),^(t)], which is a functional!) is given by 



1 



-oo 



w{X,p) = - / e'^^/^p{X + e/2,X - e/2)rfe, (4.5) 

TT J-oo 

where we see that X and ^ are our usual variables, and p has units of momentum. w{X, p) is 
somewhat analogous to a probability distribution on phase space. It is real, and integrates 
to a total of 1; its primary difference from a classical phase space distribution is that it is 
not strictly non-negative. 

If we want to advance w{X,p) in time, we can define a transfer matrix T such that 
w{tf) = Tw{to). More explicitly this is 

w{Xf,Pf) = J dXo J dpo T{Xf,pf;Xo,po)w{Xo,po) (4.6) 

where T(X/,p/; Xo,po) is defined 

T{Xf,pf-Xo,po) = di,difS^f'Pi-^<^'P^^'''f{Xf + if/2,Xf-if/2-X,^ 

(4.7a) 

f(x/,xV;xo,a;'o) = j 5x5a;' exp ^|5sys[x(t)] - 5,ys[x'(t)] + W^[a;(t), x'(t)]| (4.7b) 
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If we let t f — to = 2tt/ujo then time advancement can be performed by repeated applications 
of T. This is a sort of quantum map, 



Wi 



(4.8) 



We can ask if repeated apphcations of T will tend to converge to some invariant Wigner 
distribution Wi^v = Twi^v, analogous to the invariant measure Pmv of section III. Preliminary 



numerical calculations seem to show that this is the case |TT[. This winv appears unique, 
and should be analytic, thanks to the "blurring" effect of quantum noise. 

To make closer contact with the classical system, we might wish instead to consider 
histories in which a trajectory passes through small cells in phase space, rather than just 
intervals in X. We can write such a history by considering projections onto intervals in 
X followed very briefly by projections onto intervals in p. Histories of this type have been 
considered by Gell-Mann and Hartle, and by Halliwell, who wrote down an explicit equation 
for such a history [^|12 . 



One cannot in general specify both the momentum and position of a particle at the same 
instant. One can, however, consider a measurement of position followed by a measurement 
of momentum, and let the time between them go to zero. Halliwell calculated the probability 
of such a history using approximate projections 



1 



cixexp 



(X — X) 



|x)(a;| 



(4.9a) 



dpexp 



{p - pY 



\p) {p\. 



(4.9b) 



These are approximate projections into intervals of width and Ap of width ap, centered 
on X and p, respectively. Halliwell shows |T2| that for an initial Wigner distribution w{X,p), 
the probability of finding a particle in the phase-space cell delimited by the two above 
projections is 



p{A^,Ap 



dXdp w{X,p) exp[—a{p — p) — b{X — x 



(4.10) 



where 



when the X projection preceeds the p projection, and 

2 a? 

when the p projection preceeds the X projection. There is a restriction on these projections 
that 

< a6 < ^ ^ olal > h^. (4.13) 

To calculate the probabilities of an orbit passing through a series of such cells (Aj, Aj) 
at times tj we make use of the transition matrix T. Let us assume that the X projection 
comes first. Then it turns out that 



p({A,„ A,}) = I d{X}d{p}w{Xo,Po) 



xT(X2,p2;Xi,Pi) 

X exp[- A(p3 _ p^y -^_{X,- x,f -^{P2- PS? - ^(^2 - X.f] 

xr(X4,P4;X3,p3) 

X--- (4.14) 

where the interval Aj is centered on Xi and Aj on pi. The final projection will be of the form 
( [4.1(j| ). The case where the p projection precedes the x is very similar to the above. 

There are other ways of considering phase space projections, using coherent states or 
Gaussian combinations of coherent states. A fuller discussion of phase space histories and 



their decoherence deserves a fuller discussion elsewhere [13 . 
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B. The Master Equation. 



Another common method of studying systems such as this is by means of a Master 
equation formahsm. Caldeira and Leggett derive such an equation in the case of a harmonic 
oscillator interacting with a thermal bath at relatively high temperature . Their result is 
readily adapted to the present case, yielding the equation 

— (x, x') = - x'f'P + - x') cos{ujot)p - -^(Vix) - V{x'))p 

+2r(x - x') u 15) 

In examining this equation, the meanings of the different terms are highly intuitive. The 
K/% term is a diffusive effect resulting from the quantum noise; the F term includes the 
effects of dissipation; the q cos(to'ot) is the driving force. 

Changing to the variables X and S,-, the Master equation becomes 

%iX<& = - WAjfV + (2!g/fi){cos(^oi)(i - |(l'(A' +«/2) - V(X - il2))p 

+ (4.16) 

If ^ is small, then we can expand the potential term to give us 

ViX + e/2) - V{X - m ^ + 1^1^ + ■ ■ ■ • (4-17) 

For the Duffing potential, of course, the higher-order terms vanish. This is a convenient 
benefit of dealing with a polynomial potential. 

Transforming this equation by ( |4.5| ) gives us a new equation for the evolution of the 
Wigner distribution itself: 

^(X,p) = HK^ - ,cos(.„*)- + -(X)- 

h'd^V^^^^dhn^a^_T^^ (4.18) 



24 dX^ dp^ dp m dX 

This is almost exactly the form of the Fokker-Planck equation for the classical equation of 
motion (|2.32| ), with the diffusive hK term representing the effects of the random fluctuations 
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on the "probability" distribution and the third-derivative term being a purely quantum- 
mechanical addition, enabling w{X,p) to become negative in limited regions of phase space. 
To interpret this distribution as a probability, we must coarse-grain by averaging it over small 
volumes of phase space, producing the sort of "smeared" Wigner distribution discussed by 
Halliwell fg. 

This is not, of course, the full story. In order to correctly describe this system, we need 
not only the time evolution of the Wigner distribution, but also to specify a set of decoherent 
histories, as discussed in the previous section. Without those histories, it is impossible to 
assign classical probabilities in a consistent manner. These two approaches can be made to 
complement each other, however, as the Master equation can be solved to yield the transfer 
matrix T, defined in the previous section as a path integral. In the case of chaos, one can 
in general only solve these equations numerically, and the Master equation formalism then 
has a compuational advantage over the path integral form. 

V. INTERPRETATION OF CLASSICAL QUANTITIES 

From the previous section, we see that the behavior of a system such as we are examining 
can be evaluated on many levels: 

1. The Classical level. In the previous section we saw that all histories which deviate 
too far from the classical solution have their probabilities highly suppressed. If a system is 
large enough in scale, with enough inertia that the quantum effects are lost in other sources 
of uncertainty, we can treat it as approximately classical. Clearly, in a chaotic system this 
quantum noise does cause large alterations in the overall behavior of the system, but it is 
often impossible to separate this from thermal noise or other sources of error. 

2. The Quasiclassical level. Here we again treat the system as essentially classical, but 
now explicitly include the noise arising from quantum effects, which is large enough to be 
noticed on the scale of the system; this is the system as described by equation ( |2.32| ). From a 
practical point of view, this is the level at which quantum effects are most easily calculated. 
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This also overlaps the considerable work that has been done on dynamical systems with 
noise 0,^^-|T6[. Note that this is not the same as a semiclassical approximation, such as 
WKB. 

3. The Quantum level. Fundamentally, we can consider the system in terms of coarse- 
grainings and decoherent histories. Instead of treating a system as basically classical with 
added noise, we consider all possible histories, and compute expectation values for classical 
quantities from the probabihties of those histories. 

We've already discussed the classical (level 1) definitions of the lyapunov exponents and 
fractal dimensions used to characterize chaotic systems and strange attractors in section 
III. As pointed out, these quantities are usually defined at least formally by calculating a 
quantity for the system at different levels of coarse-graining (i.e., different box sizes e), and 
taking the limit as we go to finer and finer scales. While this has great mathematical power 
and consistency, in actual physical systems it inevitably breaks down. As Benoit Mandelbrot 
wrote on the problem of measuring coastlines with seemingly infinite levels of detail ||17|| , 

"To obtain a [fractal] Koch curve, the cascade of smaller and smaller new promontories 
is pushed to infinity, but in Nature every cascade must stop or change character. While 
endless promontories may exist, the notion that they are self-similar can apply only between 
certain limits. Below the lower limit, the concept of coastline ceases to belong to geography. 

"It is therefore reasonable to view the real coastline as involving two cutoff scales. Its 
outer cutoff Q might be the diameter of the smallest circle encompassing an island, or 
perhaps a continent, and the inner cutoff e might be the twenty meters mentioned. . . Actual 
numerical values are hard to pinpoint, but the need for cutoffs is unquestionable." 

As we shall see, in the case of chaotic strange attractors, the underlying quantum physics 
effectively provides that lower cutoff. 
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A. Lyapunov Exponents. 



Classically, the Lyapunov exponents characterize the rate at which nearby trajectories 
diverge as they evolve according to the equations of motion. In a chaotic system, one 
expects any two trajectories, no matter how close they start, to eventually move on the 
strange attractor completely independently of each other. This is measured in the classical 
case by taking the limit as points start arbitrarily near each other and evolve for arbitrarily 
long lengths of time. 

When we allow for the presence of quantum effects, however, this definition is no longer 
meaningful. As points begin closer and closer to each other, the effects of noise become 
larger and larger; one would expect the largest exponent to diverge in the limit as e — > 0. As 
we saw in section IV, the phase space cells in a decoherent history cannot be smaller than 
a certain size. This limit provides the lower cutoff mentioned above. 

For most systems it is impossible to calculate the values of Lyapunov exponents exactly. 
Instead, one performs a numerical calculation. It is easiest to calculate the highest exponent; 
lower exponents are more difficult, as their effects tend to be swamped by Ai. In a numerical 
calculation small errors are unavoidable; each such error will add a small admixture of the 
most rapidly growing component, which will quickly drown out other effects. 

Because of this, we'll first consider only the value of Ai. A simple way of estimating Ai 
is to numerically integrate equation (|2.32|) for a longer period of time, to generate a large 
number of points {xi,pi} in phase space. One can then locate nearby points, closer than a 
certain cutoff e, and trace their trajectories until they diverge further than an upper cutoff 
A. One then calculates the logarithm of the average divergence rate and averages it over 
many such pairs of points. 

I have calculated this quantity in the quasiclassical case (see fig. 4). It turns out that the 
result one calculates is not very sensitive to the upper cutoff A, but is highly sensitive to 
the lower cutoff e. In figure 4 we see the measured value of Ai as a function of e for several 
different relative strengths of the quantum noise hK. Notice how for hK > Ai diverges as 
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Because of dissipation, the overall phase-space volume of an initial distribution tends 
to decrease with time. This indicates that, classically, Ai + A2 < 0. At very small length 
scales, however, the effects of noise counteract the effects of dissipation, causing phase-space 
volume to grow rather than shrink. Thus, at small length scales we expect to see the sum 
Ai -|- A2 become positive, and eventually approach A1/A2 ~ 1; the dimension at that length 
scales should also approach an integer (2 in this case). 

We can try to define a quantum-mechanical analog of Ai. While I am not sure exactly 
what form such a definition should take, I can make a conjecture. Suppose that we start 
from the invariant Wigner distribution Win^ at t = to- We divide phase space into small cells 
{cj}, centered on average positions {vi} in phase space. These cells have a characteristic size 
e (or area e^). Let dij be the distance between the centers of the ith and jth cells. We can 
define pi to be the probability that the system is in cell q at time ti, using equation ( [4. 101) , 



and Pij to be the probability that the system is in Cj at time ti and Cj at time t2, as shown 



in ( |4.14| ). Clearly 

i j 

The probability of the system being in Cj at time t2 given that it was in q at ti is 

pm = ^- (5.1) 

Pi 

A rough estimate of the rate of spreading is then given by 

(5.2) 



A(e)qm = Pi log 



2 S PijPik(Pjk 



ph'^ 



It isn't clear whether this will agree with the usual definition of Ai in the limit. This is 
much more a rate of expansion averaged over the attractor, whereas Ai is usually defined 
as the rate at which nearby solutions diverge when followed for a long period of time. 
This latter definition has serious problems in the quantum case, where it is impossible 
to start solutions arbitrarily close together, and hence equally impossible to follow them 
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for arbitrarily long periods of time without global processes (such as folding) becoming 
important. When quantum effects are very small one can approach this long-orbit definition, 
but in that case one is really doing a quasiclassical calculation (like the one above in figure 
4), where the system can be treated as a classical stochastic equation. 

Numerical experiments might serve to explore the connections, if any, between these 
classical and quantum ideas of Lyapunov exponents. I hope to do more such exploration 
soon. Also, it is not clear to me exactly what form quantum equivalents to lower Lyapunov 
exponents might take, nor even if such a concept is useful. These questions will soon, I hope, 
have at least tentative answers. 



The information dimension is, as we saw in section IV, another number used to charac- 
terize strange attractors. It is usually defined by a box-counting algorithm of the type given 



the size of a box to go to zero no longer makes much sense. 

Instead, let us consider the information dimension Di as a function of box-size: 

EPi In Pi 



The probability pi is defined as before. For non-zero e, the exact value of Dj{e) will vary 
slightly depending on how the boxes are chosen. This ambiguity can be eliminated by taking 
Dj{e) to be the minimum possible value over all possible arrangements of boxes. In practice, 
this makes little difference. The usual classical limit is then just the limit of Dj{e) as e — > 0. 

Figure 5 shows the calculated values of Di{e) for different values of hK. As we see, at 
large length scales the fractal nature of the attractor is not readily apparent; as we shrink 
our scale, the dimension decreases, until when dropping below the lower cutoff given by the 
quantum effects it abruptly turns upward again. This was calculated quasiclassically, using 
a long orbit with associated noise. 



B. Information Dimension. 




when one takes quantum mechanics into account, allowing 




(5.3) 
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One can also do an analogous calculation using the complete quantum theory. Consider 
the invariant Wigner distribution Wi^v{X,p), defined in section IV. We can define the in- 
formation dimension -D/(e)qni using the same definition ( [5.3| ). We divide phase space into 
evenly-sized cells {q} of size e, just as in the discussion of A(e)qm above, and use the expres- 
sion ( [4. 101 ) for the probability Pi of being in the cell Cj. Again, we can eliminate ambiguity 
by minimizing -D/(e)qm over all possible divisions into cells. Clearly such a dimension will 
not even be well defined for cells of volume less than h, and will in general depend on the 
scale of the coarse- graining, just as in the quasiclassical treatment. 

Needless to say, it is much easier to extend the information dimension Dj to a proba- 
bilistic theory than it is to find an analogy for the capacity dimension Dc- The presence of 
noise will give a small but non-zero probability of finding a point in any cell, even if it is 
far from the classical strange attractor. So, for quantum dissipative chaos at any rate, Dj 
seems to be the more useful quantity. 



VI. QUANTUM CHAOS 

Since the discovery of chaos in the 1970's, there have been numerous attempts to look for 
the existence of chaos in quantum mechanical systems. Almost all of these have concentrated 
on quantized versions of non-integrable Hamiltonian systems. A seeming paradox was at the 
heart of the debate: chaos, as a classical phenomenon, depended entirely on the existence 
of nonlinear terms in the equations of motion; yet quantum systems are by their nature 
completely linear. 

In fact, this argument is clearly invalid. While it is true that the linearity of the 
Schrodinger equation and its relativistic generalizations implies that one would not expect 
chaotic behavior in the wave function itself, this has little bearing on what one would actu- 
ally see if one observed such a system. One does not measure wave functions; one measures 
particles. 

An analogous classical treatment would be to consider probability distributions in phase 
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space rather than values of position and momentum. One can then go from a set of nonhnear 
ordinary differential equations to a Fokker-Planck partial differential equation. The P.D.E. 
is completely linear. Does this then imply that chaos cannot exist in classical mechanics? 
Such a conclusion would be absurd. 

Of course, wave functions are not probability distributions, so the comparison is a bit 
misleading. A sufficiently coarse-grained Wigner distribution, however, can be made to look 
very much like a probability distribution, and its Master equation closely resembles the 
Fokker-Planck equation, as we have seen, so comparing the two is not completely inappro- 
priate. 



In fact, certain quantum models can exhibit chaos [|T^], but they are exceptional. Most 
systems which occupy a bounded volume in phase space do eventually exhibit quantum 
recurrence in which the expectation values of quantities such as energy are almost periodic 
||T9|^2[. A simple information-theory argument can be made for why this should be: a 
bounded volume in phase space V represents a finite number of possible states N V/h, 
so one would expect the long-term behavior to be periodic or quasiperiodic. According to 
Ehrenfest's theorem, a narrow wave packet will tend to follow the classical chaotic trajectory 
for a time until it has spread out to a size comparable to the phase space volume. Thus, 
it is argued, one should observe a long chaotic "transient," ending ultimately in periodic 
or almost periodic motion. As the system becomes more "macroscopic" the phase space 
volume becomes larger with respect to h and becomes longer. 



Chirikov, Izrailev and Shepelyansky describe the usual approach to quantum chaos pS 
They separate the problem into two parts, the dynamics of the undisturbed wave function 
and the effects of measurement and wave function collapse. Obviously, this periodic behavior 
of expectation values says nothing about what an experimenter would actually observe upon 
measuring the system. They compare the former with deterministic behavior and the latter 
with randomness and "noise". Clearly, an actual series of measurements would not be 
periodic at all, but would instead resemble a random chaotic trajectory. They further 
dismiss the study of dissipative chaos as a mere phenomenological approximation to an 
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underlying Hamiltonian system, e.g., in our case including both the system and reservoir 
degrees of freedom. 

If the work of Gell-Mann, Hartle, and others is to be believed, we should consider only 
decoherent histories. Since Hamiltonian systems almost by definition do not interact with 
outside degrees of freedom, one cannot really talk of "measuring" them. If measurements 
are taken, the system is disturbed; if measurements are not taken, the system evolves undis- 
turbed, but detailed histories of the motion will not decohere. In dissipative systems, by 
contrast, the chaotic system is interacting continually with the neglected degrees of freedom 
of the reservoir. These serve to provide a continual "measurement" of the system, in addi- 
tion to causing dissipation and noise, so that histories of the system variables do decohere, 
as we have seen. Thus, considering these dissipative systems from a decoherence functional 
point of view is entirely appropriate. Also, any description of real macroscopic systems must 
allow for coarse-graining over the many degrees of freedom of which we are ignorant. As one 
goes from the classical to the quantum realms, deterministic uncertainty or randomness is 
replaced by probabilistic uncertainty. 

In the Duffing oscillator model studied here, an arbitary initial distribution wq will tend to 
converge onto the invariant distribution, lim T^Wq Winv W^- This invariant distribution 



is periodic, with period 27!- /uq, so that Twinv = Winv This resembles the idea of a long 
chaotic "transient" leading to a non-chaotic, periodic behavior. 

This is not, however, very different from the behavior of classical systems evolving ac- 
cording to the Fokker-Planck equation. Most initial distributions Pq of non-zero width evolve 
into the invariant probability measure Pi^^ in much the same way that w evolves into Wi^^ in 
the quantum case (see fig. 6). The classical case is complicated by the non-analyticity and 
non- uniqueness of Pmv, so that there can be a finite probability of sitting on top of some 
unstable equilibrium or periodic point. An initial distribution with non-zero width will 
never actually become the invariant measure, of course, always itself remaining an analytic 
function if it begins as one, but the difference rapidly becomes too small to measure. These 
difficulties do not exist in the quantum case, due to the presence of noise. This difference, 
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plus the necessity for coarse-graining (which requires an initial distribution to have non-zero 
width, thanks to the uncertainty principle), is what distinguishes the classical and quantum 
cases, not the periodicity of the solution per se. The important idea is not the behavior of 
Winv, but rather the probabilities of different possible histories of the system as described in 

(EH)- 

Other authors have concentrated on the rate of expansion of uncertainty in quantum 
chaotic systems, and on how chaos leads to a form of "dissipation" in the quantum wave 
function itself ||23|-^ . 



A very interesting suggestion arose in work by Weinberg on possible nonlinear generaliza- 



tions of quantum mechanics p6| , p7| . With true nonlinearity, chaos on the level of the wave 
function could exist. Weinberg points out, however, that any nonlinearity in the theory must 
be very small to be consistent with experiment, and chaos is rare in such nearly-integrable 
systems. While such a generalization of quantum theory is certainly not ruled out be exper- 



iment, it poses a number of troubling problems. As pointed out by Polchinski either 
faster-than-light transmission of information via an EPR-type experiment becomes possi- 
ble, or different branches of histories can continue to interact with each other, effectively 
making decoherence impossible. Both of these, while certainly not impossible, have abso- 
lutely no basis in experimental evidence; indeed, the latter would imply that even classical 
probabilities could undergo a sort of "interference" with each other. 

A few papers have been published on quantum dissipative chaos p9|j30[] , mostly by in- 
cluding a phenomenological term in the Schrodinger equation to model dissipation. While 
there is a fair bit of other work on the phenomenon of quantum dissipation itself |[7|j3T|-|33 



little of it has been applied to chaotic systems, as far as I am aware, and relatively little has 
been from the decoherence point of view. 

I cannot do a broad survey of an extremely active field here. A number of good books 
now exist on quantum Hamiltonian chaos p4|,|35|. 
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VII. CONCLUSIONS 



In this paper, I have only treated a single model in depth: the forced, damped Duffing 
oscillator. It is clear, though, that the techniques used are easily applied to any nonlinear 
oscillator problem. I believe that, in general, the formalism of Gell-Mann and Hartle provides 
a rigorous method for treating any classical chaotic system. 

For some systems, of course, such a treatment may be inappropriate. Dissipative chaos 
was ffist discovered in attempts to model fluid dynamics; while in principle such problems 
could be treated in this way, any quantum effects are likely to be so small as to be irrelevant. 
For nonlinear oscillators, though, the application is quite reasonable, and it is possible that 
experiments could be done in quantum optical or electronic systems which would correspond 
to classically chaotic systems of this sort. 

What is more, using this formalism provides a link between a quantum system and its 
"classical limit" which lets one define the idea of quantum chaos in a rigorous way. And 
such systems can be treated not only semi-classically, but with the full quantum laws as 
well. 

Equally important, in this quasiclassical treatment we can see how chaotic systems can 
serve to amplify the effects of small quantum fluctuations. Sensitive dependence on initial 
conditions - the hallmark of chaos - is an important idea in understanding measurement 
situations, in which minor quantum effects can become correlated with a change in macro- 
scopic variables. In complex systems with many degrees of freedom, sensitive dependence 
on initial conditions is probably the rule rather than the exception. 

A great deal remains to be done. This formalism can readily be applied to a number of 
other systems besides the Duffing oscillator. Also, it would be valuable to develop numerical 
programs for solving the Master equation, and calculating the quantities characterizing the 
attractor in the full quantum system as well as the quasiclassical limit. Some work on this 



has already been done |TT|, but it is still in a rather crude state. But the basic outlines of 
the theory are clear. Using the theory of Gell-Mann and Hartle, a rigorous treatment of 



30 



quantum dissipative chaos is finally possible. 
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Evolution of Damped, Driven Duffing Oscillator 
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Evolution of Damped, Driven Duffing Oscillator 

Fifi. 6b t = 2 Pi 
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Evolution of Damped, Driven Duffing Oscillator 

Fig. 6c t = 4 Pi 
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Evolution of Damped, Driven Duffing Oscillator 

Fifi. 6d t = 20 Pi 
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